Abstract. New near-IR long slit spectroscopic data ob- 
tained with ISAAC on VLT/ANTU (ESO/Paranal) com- 
plement and extend our previously published near-IR data 
(AUoin et al. 2001) to produce Br-^ and H2 emission line 
maps and line profile grids of the central 4"x4" region 
surrounding the central engine of NGC 1068. The see- 
ing quality together with the use of an 0.3" wide slit 
and 0.3" slit position offsets allow one to perform 2D- 
spectroscopy at a spatial resolution w 0.5". Slit orienta- 
. . tions (PA=102° and PA=12°) were chosen so as to match 
. respectively the equatorial plane and the axis of the sus- 
I pected molecular /dusty torus in NGC 1068. The selected 
,: wavelength range from 2.1 to 2.2 fj,m is suitable to detect 
and analyze the Br.^ and H2 emission lines at a spectral 
resolution corresponding to 35 km s~^. An asymmetric dis- 
tribution of H2 emission around the continuum peak is ob- 
Qv^ : served. No H2 emission is detected at the location of the 
fSJ strong 2.2 fim continuum core (coincident within error- 
,: bars with the central engine location), while two conspic- 
uous knots of H2 emission are detected at about 1" on each 
side of the central engine along PA=90° , with a projected 
velocity difference of 140kms~^: this velocity jump has 
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, been interpreted in AUoin et al. (2001) as the signature 
of a rotating disk of molecular material. From this new 
, data set, we find that only very low intensity Br-, emis- 
sion is detected at the location of the two main knots of 
H2 emission. Another knot with both H2 and Br-y cmis- 
■ sion is detected to the North of the central engine, close to 
the radio source C where the small scale radio jet is redi- 
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Q ■ rected and close to the brightest [OIII] cloud NLR-B. It 
^ , has a counterpart to the South, placed almost symmetri- 
^ ■ cally with respect to the central engine, although mainly 
• • ,: visible in the Br^ emission. The northern and southern 
. ^ ' knots appear to be related to the ionization cone. At the 
,: achieved spectral resolution, the II2 emission line profiles 
■ appear highly asymmetric with their low velocity wing 
being systematically more extended than their high veloc- 
ity wing. A simple way to account for the changes of the 
H2 line profiles (peak-shift with respect to the systemic 
velocity, width, asymmetry) over the entire 4" x4" region, 
is to consider that a radial outflow is superimposed over 
the emission of the rotating molecular disk. We present 
a model of such a kinematical configuration and compare 
our predicted II2 emission profiles to the observed ones. 
Excitation of the II2 line is briefly discussed: X-ray irra- 
diation from the central engine is found to be the most 
likely source of excitation. Given the fact that the mate- 
rial obscuring our direct view toward the central engine 
is Compton thick (Nh > lO^'^cm"^), the observed loca- 
tion of the main H2 knots at a distance of 70 pc from the 
central engine suggests that the rotating molecular disk is 
warped. 
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1. Introduction 

Over the past decade, considerable efforts has been de- 
voted to a better understanding of active galactic nu- 
clei (AGN). In particular, the use of HST (UV and visi- 
ble ranges), the use of adaptive optics (AO) on 4m-class 
telescopes (near-IR range) and the use of interferometry 
(millimeter and radio ranges) have brought a large gain 
in spatial resolution and hence many new constraints for 
AGN mod eling. Conversely, the so-called unified model 
(see Krolik 1999|) has played an important role in stimulat- 
ing the search for AGN constituents: in particular for the 
elusive pc-scale molecular/dusty torus which is thought to 
surround the central engine (black hole and accretion disk) 
and to funnel the ionizing radiation within a cone. Models 
of the infrared (IR) emission of the torus published so far, 
have explo red torus radii f rom 1 to 100 pc (Krolik & Begel- 
Pier & Krolik |1993|, Granato & Danese |1994, 



visi ble ( Capetti et al. 1997 ) & ncar-IR (Thomps on et 



al. |2001| ), fro m AO in the near- IR (M arco et al 
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Efstathiou & Rowan-Robinson 1994, Granato, Danese & 



Franceschini 1997), while the thickness, composition and 
inclination of the torus are other key-parameters in the 
modeling. Given the input energy distribution from the 
central engine, models then predict the flux spatial distri- 
bution at IR wavelengths from 2.2 /im up to 20 /im in the 
central 100 pc region of the AGN. If one chooses an AGN 
for which the spatial resolution reachable today (typically 
0.1" with the facilities mentioned above in a range from 
UV to millimetric) provides an intrinsic scale relevant to 
the torus size, a comparison between predicted IR maps 
and observed ones will be an excellent opportunity for de- 
riving stringent constraints on the model parameters. This 
calls for investigating a close-by and torus-inclined AGN: 
NGC 1068 is an ideal target. 

This galaxy is located at a distance of 14.4 Mpc, as- 
suming a Hubble constant of 75kms~^ Mpc~^. The cor- 
responding scale is 70 pc per arcsec. Many high resolution 
images have indeed been obtained over the past decade: 
decisive contributions come from the HST in the UV, 
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Tha ttc ct al. |l997| , Rouan et al. |l998| and Marco & Al- 



loin 200C ), using shift-and-add techniques or diffraction- 
limited obser vatio ns on 8m class t elesco pes in the mid-IR 
(Bock et al. 2000 , To mono et al. 2001 ), from milli meter 
interfer omet ry (Heifer & Blitz 1995| , Tacconi et al. 1997 , 
Backer 200C , Schinne rer et al. 2000 ) and radio interferom- 
etry (Gallimore et al. 1997 ). In this paper, we consider the 
compact radio source SI (nomenclature after Gallimore et 
al. 1996) as the nuclear reference, that is the central engine 
location. Several observational facts suggest the presence 
of a molecular/dusty torus in NGC 1068. The term torus 
is used here in the sense of a rotating disk-like distribu- 
tion. From inner to outer scales, let us point out the most 
relevant ones: 

— The 1 pc-size disk of ionized material seen almost edge- 
on at PA=110°, within SI, which has been detected at 
8.4 GHz with the VLBA: it is thought to tr ace the 
inner walls of the torus (Gallimore et al. 1997) 

— The core and the 15 pc-size elongated structure (at 
PA=102°), dete cted a t 2.2 /im with CFHT AO/PUEO 
(Rouan et al. 1998) and also present on 3.5 and 
4.5 / ^m im ages with ESO AO/Adonis (Marco & Al- 
loin 200C); the position of the near-IR core is found 
by Marco et al. (1997) to be coincident with SI within 
±0.05" . This conflguration is consistent with thermal 
emission from dust particles at 1500K in the core, and 
600K at a radius of 15 pc and could trace the body of 
the torus. 



— From inter ferometer maps (Backer 2000, Schinnerer 
et al. 2000) one observes, in addition to a yet unre- 
solved core of low level CO emission, two prominent 
CO emitting peaks located symmetrically 1" (70 pc) 
away from SI at PA=98° and separated in velocity 
by 100 km s^"'^: the authors interpret the CO observed 
emission in terms of a warped CO disk. 

The possible presence of a molecular/dusty torus could 
be unveiled through the analysis of the H2 roto- vibrational 
j^=l-OS(l) emission line at rest wavelength 2.122 /im. In- 
deed, excitation of molecular transitions such as H2 2- 



E. Galliano & D. AUoin; Near-IR 2D-Spectroscopy of NGC 1068 



3 



1S(1) and H2 1-0 S(l) can occur from irradiation by UV 
photons or X-rays, as well as from shocks, all ingredients 
which are found in profusion in an AGN. In addition, as 
the H2 emission arises from warm molecular gas, the study 
of the H2 lines provides physical information complemen- 
tary to that derived from CO molecular transitions. 

Molecular hydrogen was detected for the first time in 
an extragalactic source in the central region of NGC 1068 
by Thompson, Lebofsky & Rieke ( 1978| ). They also de- 
tected Br^ emission in this region. Further H2 line obser- 
vations were carried out by Hall et al. (1981), and Oliva 
& Moorwood (1990). The first attempt at resolving spa- 
tially the H2 emitting region was made by Rotaciuc et 
al. (1991): from their best-seeing map they found that 



Table 1. Summary of the observations 



the H2 line emission was very weak at the location of the 
central engine and was mainly arising from two knots of 
unequal intensity located at ~1.3" from the nucleus along 
PA=70°. A second attempt to get a < 1" resolut ion im age 
in the H2 line emission was made by Blietz et al. (|l994|) . In 
the latter study, in contrary to Rotaciuc et al. (1991), the 
H2 emission was found to arise in three knots, the bright- 
est one being coincident (within an error bar of ±0.3") 
with the strong near-IR continuum core (central engine 
location) . 

In this study, we analyze the H2 and Br^ emission lines, 
at respective rest wavelengths 2.122 and 2.166 /^m using 
ISAAC on VLT/ANTU, in its short wavelength medium 
resolution mode. Preliminary results about the H2 emis- 
sion within a small region (1.5" x 1.5") around the cen- 
tra l engi ne were already published in a Letter (AUoin et 
al. 



2001) 



In the current paper we extend the size of the 
region explored to 4" x 4" around the nucleus and provide 
H2 line profiles and H2 line emission maps over this region. 
We also analyze, similarly to H2, the Br^ line emission 
map and line profiles. From the H2 line profile analysis, 
we build a kinematical model of the warm molecular com- 
ponent in the AGN of NGC 1068. 

We present in Section 2 the relevant information for 
data collection and data processing. In Section 3 we pro- 
vide complete sets of the H2 and Br.y line profiles across 
the entire region explored. From these data sets, emission 
maps in the II2 and Br-y lines have been reconstructed 
and are compared to the map in the near-IR continuum. 
Fluxes and velocities are also given in Section 3. From the 
II2 line peak velocities and profiles, we discuss in Section 4 
the kinematics of the molecular material and we present in 
Section 5 a simple kinematical model, aiming at reproduc- 
ing the profile shapes. Possible sources of excitation of the 
II2 line in the close environment of the AGN in NGC 1068 
are discussed in Section 5. Our concluding remarks are 
presented in Section 6. 



date 


PA 


positions 


average seeing(*) 


99/08/01 


102° 


0.6"S,0.3"S,0.0", 


0.5" 






0.3"N,0.6"N 


0.5" 


99/09/17 


12° 


0.0",0.3"W,0.6"W 


0.45" 


00/12/13 


12° 


0.3"E,0.6"E 


0.75" 


00/12/14 


102° 


0.9"S,1.2"S,1.5"S 


0.9" 


01/01/06 


12° 


0.9"E,1.2"E 


0.37" 


01/02/07 


12° 


1.5"E 


0.55" 



(*) Seeing values provided in this table have been taken from 
the ISAAC instrumental log. We have independently checked 
these values on the processed images wherever the slit was 
crossing through the nuclear unresolved K continuum source. 



2.1. Instrumental setup 

Observations were performed using the SWSl short wave- 
length arm of ISAAC, in medium resolution mode (MR), 
attached to VLT/ANTU. Due to the request for very 
good seeing quality, observations were performed in ser- 
vice mode. Results discussed in this paper are based upon 
data acquired in July and August 1999, December 2000, 
and January and February 2001. A summary of the ob- 
servations is given in Table |l|. The slit width was set to 
0.3" while its length was of 2'. The slit was successively 
positioned at PA= 102° and PA=12°. The orientation 
PA=102° corresponds to the orientation of the trace of 
the to rus su ggested by the near-IR observations of Rouan 
et al. ( 199§D , where we consequently expect to detect with 
maximum amplitude the velocity signature of a rotation. 
At the perpendicular orientation, PA=12°, no rotational 
signature is to be expected. The slit positions were chosen 
to cover a spatially continuous area around the nucleus. 
Figiu'e shows the slit positions, superimposed over the 
AO M band image from Marco & AUoin ( |200(]| ). 

The spectral resolution of the MR mode is 8700 
which, at 2.15 /im, translates into a 35kms^^ resolution. 
The projected scale on the detector is 0.147" per pixel. 
Working in the IR, we used a standard nodding technique 
(408 pixels throw). The on-source individual exposure 
time is 270 seconds. Total integration times are in the 
range from 2x270 seconds to 6x270 seconds. The signal- 
to-noise ratio of the continuum reaches 100. Standard 
stars were acquired through 0.3" and 2.0" wide slits in 
order to respectively correct for sky absorption features 
and calibrate the data in flux. 



2.2. Data reduction 



2. Data: Collection and Reduction 



Successive reduction steps were applied using the IRAF 
package, except for the electrical ghost removal, for which 
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Fig. 1. Slit positions superimposed over the M band image 
from Marco & AUoin ( 2000| ). Horizontal figures are along 
PA=102° while vertical ones are along PA=12°. 



Eclipse was used. Details on IR spectra reduction and 
especially on the nodding technique (differential compar- 
ison A - B and B - A) can be found in the ISAAC data 
reduction manual, Cuby et al. ( 2000|) . The main steps are 
the following : 

(1) correction of the electrical ghost: using the Eclipse 
is_ghost routine, 

(2) simultaneous dark, bias and sky emission line sub- 
traction: following the standard nodding technique, 

(3) flat-field division: carried out after correction of the 
illumination variation along the slit, 

(4) correction of 2D distortions (slit curvature and im- 
perfect parallelism of the dispersion direction with the de- 
tector rows) using IRAF, 

(5) correction of emission line residuals (see nodding 
technique), 

(6) extraction of the spectra, 

(7) wavelength calibration using the night-sky emission 
lines. The wavelengths are taken from the night-sky spec- 
tral at las of OH emission lines in the near-IR by Rousselot 
et al. ( |2000| ), 

(8) correction of the sky absorption lines: using spectra 
of standard stars, 

(9) flux calibration performed using stel- 
lar models found in the "stellar spectra library" 



(http://www.eso.org/instrumcnts/isaac/lib/indcx.html) 



As we were interested in the spatial distribution of 
the line emission, we extracted a series of ID-spectra 
from each 2D-spectrum: 3 pixels-high i.e. (0.45" along 
the slit height) with a sliding step of 1 pixel (0.147"). 
The final spectral resolution is 2.5 A, and the wavelength 



calibration sharp enough to measure the position of 
the night-sky lines within a precision of 0.1 A. Over the 
central 0.3" xO. 3" area, the continuum is quite intense 
and produces a fringing pattern at the level of 5% 
(peak-to-peak) which cannot be fully corrected. This is 
the region where the H2 line is weak or absent. 



3. Data: Measurements and results 

3.1. The H2 and Br^ emission: intensity maps 

In each ID-spectrum we measured the total velocity- 
integrated fluxes in the H2 and Br^ lines above the con- 
tinuum. These measurements allowed us to reconstruct 
2D maps of these line emissions. Figure || compares two 
reconstructed maps in H2, obtained independently from 
the data sets at both PAs. On the left panel, the map re- 
constructed from the data set at PA=12° has been over- 
laid on top of the map reconstructed from the data set 
at PA=102°. The seeing value indicated relates to the 
PA=12° spectrum through the intense H2 knot to the 
East of the central engine (indicated by a cross) . On the 
right panel, the map reconstructed from the data set at 
PA=102° has been overlaid on top of the map recon- 
structed from the data set at PA=12°, and the indicated 
seeing is for the PA=102° spectrum through the H2 knot 
to the East. The consistency of results derived from inde- 
pendent spectral data sets and the accuracy of the map 
reconstruction technique are successfully tested through 
the similarity of the structures appearing on both maps. 
The impact of a better image quality is also obvious in the 
left panel of Figure ^ where the strongly emitting East- 
H2 knot appears more sharply, as the seeing value was then 
down to 0.37". Four particular regions are highlighted on 
Figure 2, identified and discussed in section 4.1. These 
are the East-H2 knot, the North-H2-Br-y knot, the West- 
H2 knot and the South-Br-y knot, denominated according 
to their major line contribution. 

The ISAAC H2 map can be co mpared to the early map 
presented in Rotaciuc et al. ( 1991 ), within the limitation of 
different spatial resolutions. Considering only their H2 im- 
age obtained under the best seeing conditions (Figure 2 in 
their paper), we find indeed that their results are consis- 
tent with the ISAAC reconstructed maps, both in terms of 
positioning of the H2 knots identified and in terms of their 
relative intensity. On the contrary, the H2 line map pre- 
sented by Blietz et al. ( 1994] ) is very different from that 
by Rotaciuc et al. ( |1991 ) and from the current ISAAC 
data set. While almost no H2 line is found at the loca- 
tion of the near-IR continuum core, either by Rotaciuc et 
al. ( |1991| ) or from the ISAAC data, this is the strongest 
H2 emission knot found by in Blietz et al. (1994) - in addi- 



tion to the other H2 knots which are in reasonable agree- 



ment with those detected by Blietz et al. ( 1994) and in the 
ISAAC data set. How could we explain this discrepancy? 
As the near-IR continuum core is very strong indeed, an 
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Fig. 2. Comparison between the H2 maps reconstructed using the independent data sets at PA=102° and PA=12°. 
The seeing vahie given on each figure corresponds to the seeing value while the slit was crossing the intense eastern 
H2 emitting knot. The maps have been rebinned to a smaller square-pixel to allow easier comparison between maps 
which have originally rectangular and differently oriented pixels. Similar structures appear on both maps. The impact 
of a better image quality on the sharpening of structures is also obvious. We show as well on these maps the spatial 
areas over which the total H2, and if applicable Br-y, line fluxes have been measured and given in Table 2. These are, 
starting from the box to the East, an d tur ning clock-wise: the East-H2 knot, the North-H2-Br-y knot, the West-H2 knot 
and the South-Br^ knot (see section 4.1 for a detailed discussion). The location of the central engine is shown by a 
cross. 



extremely careful continuum subtraction is mandatory if 
one uses narrow-band imaging. The observed difference in 
H2 line emission at the position of the continuum core, 
could easily be understood if the continuum had not been 
subtracted properly in the Fabry-Perot data set by Blietz 



et al. (1994), leaving a residual as a fake H2 line emission. 
Before discussing the ISAAC maps in detail, it is impor- 
tant to assess the level of the uncertainties in positioning 
and flux measurements. 



3.2. Estimation of the error-bars in positioning and flux 
measurements 

As mentioned above, we extracted from the 2D-spectra se- 
ries of ID-spectra through 0.45" (3 pixels along the slit) x 
0.3" (slit width) windows, with a sliding step of one pixel 
(0.15") along the slit. Each ID-spectrum is assigned to 
a pseudo pixel on the reconstructed image. Hence, the 
pseudo pixel size is 0.3" (slit width) xO. 15" (shding step 
along the slit), that is three times smaller than the ex- 
traction window. This is equivalent to applying a median 
filter to an image made of independent pixels. Each 2D- 
spectrum enables to reconstruct an 0.3" wide strip of the 



final reconstructed image. These independent strips have 
to be carefully positioned along the direction of the slit 
with respect to a common reference position. Along each 
slit, the reference is set to the spectrum with maximum 
near-IR continuum level. We have checked that this posi- 
tioning procedure works fine, using an acquisition image 
of the continuum in K. For cuts along either PA=102° or 
PA=12°, we find that the near-IR continuum maxima are 
located within an 0.3" box around the intersection of the 
cut direction and its perpendicular passing through the 
core. This indicates that the relative positioning of the 
strips with respect to the core is precise within better 
than 0.3". Moreover, as the telescope offsets are set with 
a precision greater than ±0.05" and as the slit is 0.3" wide, 
in the worst case, any identified structure on the recon- 
structed map has its coordinate perpendicular to the slit 
known better than within an 0.3" window. In conclusion, 
any feature on the reconstructed images can be positioned 
with respect to the core within a box 0.3" xO.3". 

Regarding flux uncertainties, let us recall that the line 
fluxes are measured from the reconstructed line maps. To 
derive these maps we have combined spectra obtained un- 
der slightly different image quality: this introduces an un- 



6 



E. Galliano & D. AUoin: Near-IR 2D-Spectroscopy of NGC 1068 



certainty in the flux measurement, especially when mea- 
suring line fluxes from small scale structures on the im- 
age or from features with steep spatial intensity gradients. 
Actually, this effect is minimized here for the following 
reason: we have measured the velocity-integrated H2 line 
flux in the East-H2 knot, where the intensity gradient is 
steep, from a set of combined spectra all obtained un- 
der extremely good seeing (0.37"). In the other identified 
knots with H2 line emission, West-H2 knot and North- 
H2-Br^ knot, the intensity gradient is not as steep and 
then the effect of slight image quality differences not im- 
portant. On the other hand, the great advantages of our 
technique over using direct images acquired through filters 
centered on the lines of interest are first, that the contin- 
uum subtraction is very precise and second, that we avoid 
any loss of line emission resulting from an unknown and 
unexpected velocity field of the emitting material, since 
the line profile is accessible. Therefore, the final precision 
on flux measurements depends on the flux calibration step 
itself which, in this wavelength domain and given the very 
good observing conditions, is expected to be better than 
10%. 

3. 3. Line flux measurements and comparison with 
previous measurements 

We provide in Table |^ the measured line fluxes in H2 and 
Br.y of the structures identified in Figure |[ East-H2 knot, 
West-H2 knot, North-H2-Br-y knot and South-Br.^ knot. 
The "apertures" used to measure the fluxes are shown 
by the boxes sketched on Figure |[ It should be noted 
that the exposure time indicated in the headers of files 
related to the ISAAC observations of July and August 
1999 are incorrect. Consequently, the fiux measurement in 
the II2 line published by AUoin et al. (2001) is higher than 
the real value by a factor 6. The correct flux values are 
those given in the present paper. This error in the headers 
of ISAAC data files has been corrected for data taken after 
March 2000, but remains so far in the headers for data 
prior to this date. This will be modified in the future. 
Considering the fiux calibration uncertainty (±10%) and 
the possible slight flux losses due to seeing effects and 
estimated to be less than 5%, we get a final flux precision 
of around ±10%. 

For comparison with previously published measure- 
ments, those are recalled in Table |[ The dispersion in 
the H2 and Br-y line fluxes measured by various groups 
is large and there appears to be a monotonic decrease 
with time. At least two reasons can be invoked to explain 
the spread: (a) aperture differences and (b) intrinsic vari- 
ability of the AGN in NGC 1068. Regarding the second 
point, it is well established that the near-IR continuum 
luminosity of the AGN in NGC 1068 has been increasing 



(North- Il2-Br.y knot): these regions are today in a state 
of excitation corresponding to the level of activity of the 
central engine about 100 years ago. Conversely, the near- 
IR continuum luminosity is known to arise from a core 
less than 0.12" (8pc) according to Rouan et al. (1998) 
or even as small as 0.03" (0.2 pc) according to speckle 
measurements by Wittkowski et al. (1998): it may have 



increased over the past 25 years, without impacting yet 
on the H2 line emission. This may explain the observed 
decoupled behaviour between the near-IR continuum fiux 
and II2 line time variations. Would well-sampled light- 
curves be available in the near-IR continuum and H2 line, 
then time-delay effects could be used indeed to probe in 
another way the size of the H2 line emitting region (e.g. 
Peterson |l99l . 

We have measured as well the continuum value at 
2.18 ^m on the ISAAC spectra passing through the core. 
Taking into account the measured seeing (FWHM of a cut 
through the continuum), considering two values for the 
intrinsic core size of respectively 0.12" and 0.2" FWHM, 
and given the 0.3" slit width, we have computed the cor- 
responding slit losses. The final values for the flux of 
the core continuum at 2.18 /zm are respectively 2.06 10^^" 
and 2.16 lO'^*^ erg s'^ cm" Vni"\ for an 0.12" and 0.2" 



FWHM core. The value derived by Rouan et al. ( |1998D 
from AO observations of the unresolved core in K was of 
0.8 evgs^^ cmr"^ HTcT^ . 



Table 2. Line fluxes for the identified knots 



Source 


Aperture 
center 


Aperture 
size 


H2 line flux 


East-H2 
West-Ha 
South-H2 
North-H2-Br^ 
Total 


±1.0", ±0.3" 
-1.3" ,-0.3" 
0.2" ,-0.4" 

-0.1" ,±0.7" 
0",0" 


1.0" X 1.4" 
1.0" X 1.4" 
0.8"x0.8" 
1.2" X 1.0" 
4.0"x4.0" 


2.3±0.3 
0.7±0.1 
0.15 ±0.05 
1.3±0.2 
6.1±1.0 


Source 


Aperture 
center 


Aperture 
size 


Br^ line flux 


East-H2 
West-Ha 
South-H2 
North-H2-BrT, 
Total 


±1.0", ±0.3" 
-1.3" ,-0.3" 
0.2" ,-0.4" 

-0.1", ±0.7" 
0",0" 


1.0" X 1.4" 
1.0" X 1.4" 
0.8"x0.8" 
1.2" X 1.0" 
4.0"x4.0" 


0.25 ±0.15 
< 0.2 
0.4 ±0.1 
2.8 ±0.5 
5.2±1.0 



steadily since 1974 (Marco & AUoin p000| . Glass |1997D . 
However, notice that the H2 line emission is arising from 
regions located at least 40 pc away from the central engine 



Note: The flux unit is 10~ ergs" cm~ . The apertures are 
shown on Figure ^ ; the aperture centers are deflned with refer- 
ence to the position of the K band unresolved core, coincident - 
within error-bars - with the radio source SI (coordinates J2000 
RA=0242 40.7098, Dec= - 00 00 47.938 after Bland-Hawthorn 
et al. 1997). Be aware that the uncertainty on this positional 



reference is around 0.1" 
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Table 3. Previously published fluxes in the H2 and 
Br^ lines 



Author 


Aperture 


Ha flux 


Br^ flux 


Hall (1981) 


3.8" diam. 


16.5 


18.5 


Oliva (1990) 


6" diam. 


15. ± 1 


13. ± 2 


Rotaciuc (1991) 


~ 3" diam. 


8.6 


4. 


this paper (2002) 


totals 3" diam. 


6.1 ± 1.0 


5.2 ±0.5 



The flux unit is 10 erg s ^ cm ^ . If not mentioned, the error- 
bars were not provided by the authors. 



3.4- The H2 and Br^ emission: line profiles 

Figure ^ and Figure ^ provide respectively grids of 
the H2 and Br..^ line profiles. In each box across this 
grid, we sample the line profile and flux from an 
0.3" X 0.45" sub-area. The X-axis provides its offset along 
PA=102° (close to East- West) and the Y-axis its offset 
along PA= 12° (close to North-South). Positive coordinates 
along both axes are respectively to the East and North. 
The unresolved continuum core in the K band is referenced 
at coordinates (0.0", 0.0"). Regarding the central -0.3" to 
-1-0.3" strip, the strong near-IR continuum from the core 
induces a residual fringing level (3%) and a noise level 
which prevent measuring weak H2 line emission. Given 
the mean seeing value (0.5" in the K band), line profiles 
in contiguous boxes are not fully independent one from 
the other. As the goal of this study is to understand and 
model the overall kinematics of the molecular gas across 
the studied region rather than its very local behaviour on 
a scale smaller than 0.5", this inter-dependence does not 
limit much the interpretation. 



coincides with the [OIII] emitting cloud NLR-B seen on 
the HST maps. 

We observe that the H2 line emission is not distributed 
uniformly around the central engine. It consists of a low 
level, extended and diffuse component which is minimum 
where the continuum is maximum (central near-IR core 
at (0.0", 0.0")). In addition to this extended H2 line emis- 
sion, three knots which are particularly intense in the 
H2 line emission have been identified: the correspond- 
ing line profiles (H2 and/or Br-^) are displayed in the in- 
serts in Figure ||. The most conspicuous one is the East- 
H2 knot, well defined in position and velocity (insert A). 
We can isolate as well the West-H2 knot (insert B), al- 
though it appears weaker in intensity (factor around 3). 
No Br-y emission is detected at the location of the West- 
H2 knot, while the East-H2 knot does show a low level of 
Br^ emission. This emission can be attributed to the wings 
of the strong Br.y emission coming from the North-H2- 
Br^ knot, rather than being intrinsic to the East-H2-knot 
itself. We notice that, within positional uncertainties, the 
East-H2 knot and the West-Il2 knot are located almost 
symmetrically with respect to the position of the central 
engine. Their H2 line peak velocities are respectively at 
1074±4kms~^ and 1214±4kms~^. Under the assumption 
they are kinematically associated in the form of a struc- 
ture rotating around the central engine, this leads to a 
rest velocity of the central engine of 1144±4kms~^. This 
velocity is fully consistent with the HI heliocentric veloc- 



4. Kinematics and distribution of the emitting H2 
material and of the ionized hydrogen over the 
central 4" X 4" area 

Emitting H2 material and ionized hydrogen over the 
central 4"^ 4" area 

Figure ^ summarizes some of the conspicuous facts to be 
concluded from these observations: superimposed over the 
map in the H2 emission line, the contours indicate the flux 
in the Br-^ emission line, while the cross features the lo- 
cation of the K band continuum maximum on our ISAAC central engine. Conversely, in the North-Il2-Br^ knot, the 

blue-peak of the H2 line profile is unique to this location, 
which is spatially coincident with the strongest Br^ line 
emission detected. Although at this location the Br-y line 
does not show an obvious double-peaked profile in cor- 
respondence to that observed in the H2 line, one cannot 
exclude that it might be the possible blend of two compo- 
nents with greater FWHM (^700kms~^) than the II2 line 



ity of the host galaxy, 1148±5 kms"^(Brinks et al. |T99j. 
In addition, the H2 line profiles in the East-H2 knot and 
West-Il2 knot have similar width and are symmetrical in 
shape relative to the galaxy rest velocity. The line profile 
of the East-Il2 knot has a blue-shifted peak and shows an 
extended red wing, i.e. toward lower velocities. The line 
profile of the West-H2 knot has a red-shifted peak and 
also displays a wing toward lower velocities (extended blue 
wing). The H2 line profiles (shape and velocity) change 
continuously and smoothly all over the region of diffuse 
emission and the identified East-H2 and West-H2 knots: 
hence we can reasonably assume that they belong to a 
same kinematical component. The third knot identified is 
located to the North (insert C). It shows different char- 
acteristics than the East-H2 knot and West-Il2 knot and 
emits both in H2 and Br^. In this North-Il2-Br-y knot, the 
H2 line profile is double-peaked. If one considers only the 
red peak of the II2 line in the North-H2-Br-), knot, one finds 
that the II2 line profile peak velocity varies in an organized 
fashion from the eastern side to the western side of the 



data set. It has been shown by Marco et al. (1997) and 
Rouan et al. ( |199^ ) from AO images in the K band, that 
this unresolved core (less than 0.12 pc in size) is coinci- 
dent, within error-bars, with the radio source SI (cen- 
tral engine). The triangle shows the location of t he ra - 
dio source C (Muxlow et al. 1996 ; Gallimore et al. 1996|) . 
Marco et al. ( 1997 ) also showed that the radio source C 
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Fig. 3. The H2 emission line profiles across the 3"x3" region around the nucleus. North is up, East is left. Each 
box corresponds to an emitting region of 0.45" x 0.3". The unresolved continuum core in the K band is located at 
(0.0" ,0.0"). The X-axis corresponds to PA=102° and the Y-axis to PA=12°. The velocity scale of the profile is shown 
on the bottom left of the X-axis. The intensity scale is displayed on the bottom left of the Y-axis. 



components (~190 kms""'^). Indeed the Br^ line emission in Br.^ than H2. Notice that the Br^ emitting knots, the 
is expected to arise from the NLR where the presence of North-H2-Br^ knot and South-Br^ knot, are also quasi- 
ionized material at high velocity is inferred from the ob- symmetrical with respect to the central engine and outline 
servation of broad wings in the [OIII] line profiles (AUoin the two sides of the ionizing cone. 

et al. 1982, Cecil et al. lOQCf ). Let us remind finally the From our results we can infer the following structure of 

South-Br-y knot which appears to be a stronger emitter the molecular and ionized material: 
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Fig. 4. The Br-y and H2 emission line profiles across the 1.2" xl.2" region to the North of the nucleus, corresponding 
to the North-H2-Br^ knot. The Br^ line profiles are shown with a red and thick line while the H2 profiles are shown 
with a black and thin line. Note that the X-axis velocity scale had to be changed because of much wider profiles in 
Br^ than in H2. As in figure |[ the X-axis corresponds to PA=102° and the Y-axis to PA=12°. The orientation is the 
same as in figure |[ 



The molecular H2 line emission from (a) the diffuse 
component, (b) the East-H2 knot, (c) the West- 
H2 knot and (d) the red peak on the H2 profile 
from the North-H2-Br-y knot, all arise from a same 
system with organized kinematical properties. This 



is demonstrated by the continuous spatial change of 
the H2 line profile peak velocity. Within this system, 
the maximum velocity shift observed between the 
East-H2 knot and the West-H2 knot has a value of 
140 km s^^, these knots being located at a distance 
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Fig. 5. Molecular and ionized gas structures in the central 4" x4" of NGC 1068, as shown on a composite image made 
from the independent maps restored from the ISAAC data sets obtained at PA=102° and PA=12°. The grey-scale 
map corresponds to the velocity-integrated H2 line emission. We have superimposed with black contours the Br-^ line 
emission which is detected mainly to the North and South of the central engine. The labels on the contours are 
the surface brightness in 10^'^ ergs^^ cm^^ sr^^. The cross shows the location of the central engine (position of the 
unresolved core seen in the K band continuum: its coincidence with the radio source SI is ascertained within the 
error-bar of AO measurements by Marco et al. 1997). The triangle indicates the location of the radio source C. In the 
inserts A, B and C, we have displayed the H2 line profiles (thick line) and the Br-y line profile (thin line) for three 
particular regions, the East-Il2 knot, the North-Il2-Br^ knot and the West-H2 knot. The profiles are displayed on a 
velocity scale to allow an easy comparison: the velocity zero point corresponds to the rest wavelength of each of the 
lines after correction for the heliocentric recession velocity of the galaxy, 1144 ±4kms^^, measured on our data as the 
mean velocity of the two sides of the suspected torus (see section 4.1). 



of ~1" (70 pc) on each side of the central engine. 
The simplest interpretation of this velocity jump is in 
terms of the signature of a rotating structure: under 
the assumption of keplerian rotation, this implies a 
central mass of 10^ M© (AUoin et al. |2001|). Modeling 



of the kinematics of this system, derived from the 
H2 line profile changes over the entire region, is 



presented in Section 4. 

— The Br^ line emission and the blue-peak of the II2 line 
emission in the North- Il2-Br-y knot, as well as the 
Br^ emission in the South-Br.y knot belong to a kine- 
matically distinct system. This system is probably as- 
sociated with the NLR material, as suggested by the lo- 
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cation of the North-H2-Br^ knot and South-Br^ knot, 
by the absence of velocity offset of their hne profile 
peaks and by the large width of the Br^ line. In addi- 
tion, let us notice that the North-H2-Br-y knot is found 
to be close to the radio source C identified on the MER- 



H2, PA=12, 0.9"E 



PA=1 2. 



LIN map (Muxlow et al. 1996). The source C coincides 
with the re-direction of the small-scale radio jet (ini- 
tial direction PA=12°), possibly induced by the impact 
on th e NLR -B cloud (for details see Bland-Hawthorn 
et al. 1997 ). The ionization resulting from the impact 
would explain the detection of the [OIII] and Br-^ lines. 

4-.2. The molecular material over the 20"x20" area 

The ISAAC long slit data set also provides information 
about the distribution of the H2 line emission over a larger 
scale. The 2D spectra, particularly along PA=12°, show 
weak H2 and Br^ emission up to almost 20" from the 
nucleus on both sides. These two lines are most clearly 
detected between 14" and 18" from the nucleus to the 
North, and between 12" to 17" from the nucleus to the 
South. The locations of these various emission spots cor- 
respond to the star forming spiral arms clearly traced by 
th e CO millimeter lines (Tacconi et al. 1997 , Schinnerer et 
al. 200C| ). The line intensity ratio Br^/H2 is about 3 to the 
North and 1 to the South. These lines have a FWHM of 
35kms^^ and are much narrower than the emission lines 
in the central 4" x 4" region. The spatial coverage of our 
data set is too scarce to extract complete 2D information 
over the central 40" x40" where the CO arms and bar de- 
velop. Still, we notice that, along PA=12° the presence of 
H2 and Br.y is conspicuous (Figure ||), while these lines 
are rather weak along PA=102°. Along PA=12°, the ve- 
locity curve does not exhibit a simple behaviour. Along 
PA=102° the velocity curve traced by H2 and Br^ is con- 
sistent with a logarithmic potential: zero on the nucleus 
and then rising up to a plateau. 



5. Kinematics: a rotating disk and an outflow 

An important constraint on the kinematics of the molecu- 
lar material around the central engine in NGC 1068 comes 
from the velocity jump of 140kms~^ between the East- 
H2 and the West-H2 knots. An ad-hoc jet configuration 
along PA-102° might produce such projected velocities; 
however, as the ionization cone axis and the radio emission 
all appear to be aligned along PA=10°, it seems rather 
natural to interpret the observed molecular emission dis- 
tribution and kinematics in the framework of a somewhat 
axisymmetrical structure located in a plane perpendicular 
to the axis defined at PA=10°. Another strong constraint 
on the kinematics comes from the asymmetry of the ob- 
served H2 line profiles. They exhibit a more extended wing 
on the side of their small velocities: ie. blue wing for the 
profiles which are redshifted with respect to the galaxy 
rest-frame velocity, and red wing for the profiles which 
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Fig. 6. Extended H2 and Br-^ emission up to 15" from 
the nucleus and with the slit passing through the East II2 
knot. The left and right figures display respectively the 
(P-V) diagram for II2 and for Br-y, along PA=12°. 

are blueshifted. This observational fact must be accounted 
for. 

5.1. Modeling the line profiles: a model with two 
kinematical components 

The hypothesis of a rotating structure for the molecular 
material is tested through a simple model. The model is 
used to produce a grid of II2 emission line profiles to be 
compared to the observed ones over the 4" x4" observed 
region. The model allows to produce profile grids for either 
of two independent kinematical components and for their 
combination. The two kinematical components considered 
are: a rotating disk and an outflow, as illustrated on Figure 

For the rotation of the disk we assume a velocity field 
induced by a logarithmic gravitational potential. Such a 
velocity field is described by two parameters: a character- 
istic radius rp and a characteristic velocity Vp. The velocity 
curve flattens at radii greater than rp, and reaches a ve- 
locity Vp. The outflow is only parameterized through the 
outflow velocity. 

The model is built in the following way: a 2D grid, 
where each element represents a small surface of the disk, 
is placed in the plane of the sky. We choose the X-axis to 
represent PA=102° and the Y-axis to represent PA=12°. 
3D vectors are assigned to each element of the grid: one 
for the position and one for the velocity of each kinemat- 
ical component. A rotation around the X-axis, to account 
for the inclination of the disk, is applied to these vectors. 
Then a v^ "image" is built for each kinematical compo- 
nent; we assign to each pixel of this image the along-the- 
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Fig. 8. Along-the-line-of-sight Vz velocity map of the rotating component (left) and of the outflow component (right) 
in the considered model 



narrow rotating component 



broad rotating component 



tating + ouflowing components 




(a) 



-100 -50 50 

velocity [km s"'] 




-400 



-aOO 200 

velocity [km s"'] 




-400 -200 200 

velocity [km s"^] 



Fig. 9. Simulated profiles for the position [1"E,0"N] on the grid. From left to right: (a) rotating component alone, 
intrinsic line FWHM=15 kms"'^ and instrumental profile not taken into account: the profile is asymmetric because of 
the mixture of different projected velocities resulting from the seeing, (b) rotating component alone, with intrinsic hne 
FWHM=150kms~"'^: the profile is symmetric because the intrinsic width is much greater than the range of projected 
velocities in a seeing sized-area. (c) rotating plus outflowing components with intrinsic line FWHM=150kms~^ for the 
rotating component and line FWHM =200kms~^ for the outflowing component: this combination allows to reproduce 
at the same time the width and asymmetry of the observed profiles, as well as the spatial evolution of the asymmetry. 



line-of-sight velocity of the grid element whose projected 
position on the sky corresponds to the pixel position. The 
Vz images for a rotating and for an outflowing compo- 
nent are shown in Figure ^. From each image, we build 
a 3D cube where the third dimension is the velocity. A 
Gaussian profile is used to represent the emission line and 
is centered at the corresponding Vz value. The intensity 
assigned to the profiles was chosen to be constant, after 
checking that this assumption does not affect our conclu- 
sions. Each plane from this data cube is then convolved 
with a 2D Gaussian representing the point spread func- 
tion (PSF) to account for seeing effects. The parameters 
for the kinematical model are the following: 



— the velocity field parameter(s) (two parameters for the 
rotation and one for the outflow) 

— the inclination of the disk 

— the intrinsic width of the emission line for each kine- 
matical component 

— the relative intensity of the emission from each of the 
kinematical components 

The H2 line profiles predicted in the two configura- 
tions - with or without outflow - are shown in Figure ^ 
We find that it is not possible to reproduce at the same 
time the general characteristics (asymmetry, width and 
peak position) of the observed H2 line profiles using 
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Fig. 7. Illustration of the two-component kinematical 
model. Arrows represent the motion of the two compo- 
nents. 

a model with only a rotating disk. On the contrary, 
the two-component model (rotating disk plus outflow) 
reproduces well all the profile characteristics and we shall 
pursue the profile fitting with such a kinematical model. 
Do we expect, from a modeling point-of-view, such an 
outflow to be present in the environment of the AGN? 
In classical AGN models, the observed conical funneling 
of the ionizing radiation fleld is explained by inflation of 
the inner part of the molecular/dusty torus. One may 
therefore anticipate that the impact of the X and UV 
photons on the edges of the inflated inner regions of 
the torus will sweep away material from its surface and 
generate an outflow. 



5.2. Results 

The simulated H2 profiles over the 3" x 3" region are plot- 
ted on top of the observed ones in the grid shown in Fig- 
The parameter values used to build the grid of 



ure 
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simulated profiles matching best the observed ones, are 
the following: 

— Vp— lOOkms"^ 

— rp= 50 pc 

— outflow velocity: 140 km s""'^ 

— inclination angle: 65° 

— FWHM of the rotating component: 150 km s""'^ 

— FWHM of the outflowing component: 200 km s""'^ 

— ratio between the line intensities emitted by the rotat- 
ing and outflowing component: 2 

The modeled profiles have been scaled to peak intensity 
in order to allow an easier comparison among the pro- 
file shapes. The agreement between the observed and 
modeled H2 line profiles is very good over the entire 
area. This validates the choice of a two-component 
kinematical model. One should note however that 



the figures derived for the model parameters are not 
constrained in a unique manner and some variation 
around the quoted figures may be expected as line 
profiles from other molecular transitions arising in the 
warm molecular gas will become accessible to observation. 

In comparison, the CO data set discussed by Schin- 



nerer et al. (2000) has been modeled in the framework of 
a warped disk. In the CO modeling, little emphasis has 
been put on fitting the CO line profiles, while in the case 
of the H2 line emission, the profiles are used as strong con- 
straints for the modeling. Let us recall in addition that low 
order CO transitions probe cold molecular material, while 
the H2 transition studied here arises from warm molecular 
matter. Therefore a detailed comparison between the two 
types of models is hampered by the use of different con- 
straints on one hand, and by reference to different temper- 
ature phases of the molecular material, on the other hand. 
We note however that the main emission knots detected in 
H2 are located at positions coincident with the CO emit- 
ting knots. It is quite conceivable that the H2 line origi- 
nates at the surface of molecular clouds exposed to intense 
radiation fields, while the CO line arises from the denser 
and cooler inner part of these clouds. In such a configu- 
ration, surface-effects (wind...) would have a great impact 
on the shape of the H2 line profile while not affecting the 
shape of the CO line profile. 

6. Discussion: source of excitation of the H2 line 

Apart from the North-Il2-Br.y knot, the lack of strong 
Br^ emission coincident with II2 emission suggests that 
most of the molecular material emitting the H2 line is 
not directly irradiated by UV photons which would also 
ionize hydrogen. Supporting a similar conclusion is the 
90° angle offset between the ionization cone axis ( along 
the South-North direction, Macchetto et al. 1994 ) and 
the predominant direction of II2 line emission (East- 
West) in NGC 1068. Moreover, it has been shown indepen- 
dently that the observed H2 line ratio 2-1 S(1)/1-0S(1) in 
NGC 1068 is consistent with thermal em ission rather than 
with UV-pumping (Oliva and Moorwood 1990| ). Therefore, 
these three facts lead us to disregard UV photons as the 
source of excitation of the II2 line in the close environment 
of the central engine in NGC 1068. 

The remaining two possibilities to excite H2 to the re- 
quested ~2000K are, (a) heating by shock waves, or (b) 
heating due to irradiation by X-rays (Mouri 1994). 
Could shocks play a dominant role in the excitation of the 
H2 emission in the particular case of NGC 1068? Following 
Draine, Roberge & Dalgarno ( |1983| ), the shock velocities 
required are of the order of Vshock ~ 10 — 30kms~^. An 
order of magnitude of the total shock luminosity is given 
by L 

shock ^ ^I^Pcloud Vshock ^ohock-, where pdoud is the 
density of the cloud, Vshock the shock v elocity and Ashock 
the total shock area (Heckman et al. 1986). A density 
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Fig. 10. Modeled profiles (dotted lines) and observed profiles (solid lines) over the 3'x3" central area. Values of the 



model parameters can be found in section 5.2 



nH2=10^cm ^, a shock velocity Vghock—'^O km s and a 
shock area Ashock=2TTrli^^^ for rc;o«d=20pc are reason- 
able values for the East-H2 knot. The total shock luminos- 
ity for these parameters is lO^'^ergs"^. About one tenth 
of this luminosity is radiated through the 1-0 S(l) H2 hue. 
The resultant predicted luminosity in this line is thus 6 or- 
ders of magnitude below the observed 3.10"^^ erg s~^ emit- 
ted in the East-H2 knot. Even if we increase generously 



riHa or Aghock, the observed flux values cannot be ac- 
counted for in this way. Therefore we can disregard shock 
heating as the source of excitation of the H2 line around 
the central engine in NGC 1068. 

Maloney ( |1997D has already discussed the X-ray irradia- 
tion of molecular clouds by the central X-ray source in 
NGC 1068. His model was computed for an intrinsic 1- 
100 keV luminosity of 10'^'' erg s^^ for the central engine 
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and an attenuating column density of 10^^ cm~^ between 
the central engine and the H2 emitting molecular clouds. 
With this set of parameters, his model accounts for the 
observed H2 intensity and predicts the right size-scale of 
the H2 line emission. However, Matt et al. (1997) have 



shown that the opaque material which obscures our direct 
view to the central engine of NGC 1068 is Compton thick 
(Nh > lO^'^cm"^). If this opaque material has a more or 
less toroidal distribution - as suggested by the existence 
of the ionizing cone -, then the column density between 
the central engine and the molecular material found at the 
location of the East and West-H2 knots must also be of 
the order of at least lO^^'cm"^. In such a configuration, 
the model predicts an H2 emission arising much closer to 
the nucleus than the observed 70 pc. An alternative pos- 
sibility is that the azimuthal density distribution in the 
torus is such that the opening angle for the X-rays is larger 
than for the UV radiation. This situation is likely to occur 
since X-rays can penetrate much larger column densities 
than UV photons, implying that the viewing angle range 
over which the torus is Compton thick is smaller than the 
viewing angle range over which it is UV-opaque. Still, in 
such a situation, the H2 emitting clouds must be located 
within aperture of the cone inside which the X-rays can 
escape. This suggests that the molecular clouds located 
at a distance of 70 pc (where the H2 emission is mainly 
detected) are not strictly orbiting in the plane perpendic- 
ular to the axis of the opaque torus. Consequently, either 
the H2 emitting material is a component distinct from the 
torus and is simply orbiting in the plane of the galaxy 
(inclination ~ 40°), or we are dealing with a single disk 
structure within which the inclination of the plane, where 
the molecular material orbits, changes continuously with 
radius: in other words the disk is warped. 
The presence of a warped structure in the AGN of 



NGC 1068 has also been proposed by Backer (gOOOj) to fit 
both the intensities and the kinematics of the CO emitting 
material. Inter estingly, as shown in a theoretical work by 
Quillen ( 2001 ), a wind blowing across a dense, thin and 
initially flat rotating disk would automatically generate 
a warp. This suggests that there is a possible physical 
connexion between the warp and the outflow, hence be- 
tween the observed distribution of II2 emission and the 
H2 line profile shapes. First, the asymmetry of the pro- 
files witnesses the presence of an outflow. Second, model- 
ing predictions tell us that this wind would have naturally 
induced a warp of the molecular disk. Third, it is the ex- 
istence of the warp which makes possible the irradiation 
of the molecular material by the central X-rays. In sum- 
mary, this scenario makes the wind the indirect cause of 
the H2 emission being observed at a rather large distance 
from the central engine. 



7. Concluding remarks 

We have provided in this paper a comprehensive study of 
the kinematics of the warm molecular material in the en- 
vironment of the AGN in NGC 1068. New runs of 2D spec- 
troscopic observations (ISAAC on VLT/ANTU), obtained 
under very good seeing conditions, have complemented 
the data set discussed earlier by AUoin et al. ( ^001 ). 
The achieved spectral (35kms~^) and spatial (0.5") res- 
olutions allow to recover emission maps in the H2 and 
Br^ emission lines, at respective rest wavelengths 2.122 
and 2.166 /im , over the central 4"x4" region surround- 
ing the central engine in NGC 1068. In addition, one wit- 
nesses the spatial evolution of the II2 and Br^ line profiles 
over the inner 3" x 3" region with a spatial sampling of 
0.3"x0.3". 

The H2 and Br^ line emission maps and profile spatial evo- 
lution have been used to analyze the distribution and kine- 
matics of the warm molecular gas across the 200 pc central 
region. The salient results of this work are the following: 

— There is no detectable H2 emission at the location of 
the strong 2.2 /im continuum core (which is known al- 
ready to be almost coincident with the central engine) , 
while two main regions of H2 emission, East-H2 and 
West-H2, are detected at about 70 pc on each side of 
the central engine along PA=90°, with a projected ve- 
locity difference of 140kms~^: this velocity jump is 
interpreted at first order as the signature of a rotat- 
ing structure. Under this assumption, the systemic ve- 
locity of the AGN is found to be 1144±4kms"^ No 
strong Br^ emission is detected at the location of these 
knots. The East-H2 knot is about three times brighter 
in H2 line emission than the West-H2 knot 

— One detects as well a diffuse and extended H2 emission, 
although at low intensity level. At about 0.3" to the 
North of the central engine, the H2 emission also co- 
incides with a Br-y line emitting knot. This North-H2- 
Br-|, knot is close to the radio knot C and the brightest 
[OIII] cloud, NLR-B, although not exactly spatially 
coincident with them. 

— To the South, another Br^ line emitting region is 
observed, although weaker in intensity, the South- 
Br-y knot. The North-H2-Br-y knot and South-Br^ knot 
are almost symmetrical with respect to the central en- 
gine and appear to be located on the inner edges of 
the ionizing cone. 

— The observed changes in the H2 line profile across the 
200 pc central region cannot be reproduced by the sim- 
ple, disk-like, kinematical structure which matches the 
velocity jump between of the East-H2 knot and the 
West-H2 knot. A possible solution is to add an outflow 
at the surface of the rotating disk-like structure. The 
main parameters of the complete kinematical model 
are: a rotational velocity around lOOkms"^ at a char- 
acteristic radius of 50 pc , an inclination angle of the 
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disk-like structure of 65° and an outflow velocity of 
140 km s-^ 

— The source of excitation of the H2 line is briefly dis- 
cussed. Both UV photons and shocks can be discarded 
as the main source of excitation of the H2 line, while 
X-ray irradiation from the central engine is found to 
be the most likely mechanism. 

— An additional consequence of an outflow at the surface 
of a rotating disk-like structure is the occurrence of a 
warp. This warp does naturally explain how the X-rays 
can impact the molecular disk at a distance as large 
as 70 pc and give rise to the two observed main knots 
emitting H2. 
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